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Abstract. The canonical problem of the propagation of surface-plasmon-polariton 
(SPP) waves localized to the planar interface of a metal and a sculptured ne- 
matic thin film (SNTP) that is periodically nonhomogeneous along the direction 
normal to the interface was formulated. Solution of the dispersion equation ob- 
tained thereby confirmed the possibility of exciting multiple SPP waves of the 
same frequency or color. However, these SPP waves differ in phase speed, field 
structure, and the e-folding distance along the direction of propagation. 
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1 Introduction 



Among the various forms of electromagnetic surface waves, the surface-plasmon- 
polariton (SPP) wave has the longest history of theoretical development and ap- 
plication. Arising from earlier developments, the SPP wave was envisioned in the 
mid-twentieth century as a wave guided by the planar interface of a metal and an 
isotropic dielectric material. Since that time, the theory has evolved to encom- 
pass interfaces between a metal and various dielectric materials of greater com- 
plexity. The inclusion of anisotropic, homogeneous dielectrics[l, 2, 3, 4, 5, 6, 7, 8] 
in the study of electromagnetic surface waves has been considered for some time 
now. More recently, SPP waves at the interface of a metal and a nonhomoge- 
ncous dielectric material, which exhibit several interesting properties, has been 
receiving considerable attention. The nonhomogeneous dielectric materials in- 
vestigated include continuously varying materials[ll, 12, 9, 10] such as liquid 
crystals, as well as layered structures [13, 14, 15, 16]. 

This series of papers is devoted to theoretical and experimental investiga- 
tions on the propagation of a surface electromagnetic wave guided by the planar 
interface of a metal and a sculptured nematic thin film (SNTP), the latter be- 
ing periodically nonhomogeneous in the direction normal to the interface. In 
Parts I [17] and II [18], the absorbance, reflectance and transmittance of an 
incident linearly polarized plane wave were calculated when this planar inter- 
face is implemented via a metal-SNTP bilayer and a prism in a Kretschmann 
setup. The wavevector of the incident plane wave was supposed to lie wholly 
in the morphologically significant plane of the SNTP, but that restriction was 
lifted in Part II. The computed results showed that multiple surface-plasmon- 
polariton (SPP) waves, all of the same frequency or color, can be excited at the 
metal/SNTP interface. The guided SPP waves possess different field structures 
as well as different phase speeds. 

In Part III [19], experimentally obtained data was presented in support of 
the theoretical predictions of Part I. The absorbances of two different metal- 
SNTP bilayers were evaluated as the difference between unity and the measured 
reflectance, the transmittance being assumed to be null- valued at angles of inci- 
dence exceeding the critical angle for the interface of the SNTP with the prism. 
Analysis of the collected data conflrmed the possibility of exciting multiple SPP 
waves with different phase speeds and held structures. Parenthetically, in a par- 
allel effort the same possibility was theoretically predicted [9, 10] and validated 
experimentally [20] for the planar interface of a metal and a chiral sculptured 
thin fllm [21, 22]. 

The method adopted for theoretical predictions and experimental verifica- 
tion in Parts I-III was an indirect one, because the existence of SPP waves was 
deduced from the absorbance, refiectance and transmittance characteristics of 
a metal-SNTP bilayer. Each layer in the bilayer is of finite thickness. In the 
canonical problem of SPP-wave propagation, the metal and the SNTP have to 
be semi-infinite. The solution of the canonical problem provides incontrovertible 
proof of the existence of multiple SPP-wave modes, and also eliminates possible 
confusion with waveguide modes spread over the entirety of the SNTP. Accord- 
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ingly, for this paper, we set out to prove the existence of multiple SPP waves 
directly. The approach adopted by us is similar to that of Agarwal et al. [23], 
being independent of any incident plane wave. The formulated boundary-value 
problem was solved numerically for the SPP wavenumbers. 

The canonical problem is formulated in Sec. 2, and numerical results are 
presented and discussed in Sec. 3. An exp(—iiot) time-dcpendcncc is implicit, 
with uj denoting the angular frequency. The free-space wavenumbcr, the free- 
space wavelength, and the intrinsic impedance of free space are denoted by 
kf, = w^eo/io, Ao = 2tt /kg, and rja = ■\/ Mo/fo, respectively, with /io and Eq being 
the permeability and permittivity of free space. Vectors are in boldface, dyadics 
are underlined twice, column vectors are in boldface and enclosed within square 
brackets, and matrixes are underlined twice and square-bracketed. The asterisk 
denotes the complex conjugate, and the Cartesian unit vectors are identified as 
iix, Uy, and u^. 



2 Theory 

Let the half-space z < be occupied by an isotropic and homogeneous metal 
with complex- valued relative permittivity scalar Cmet- The region z > is 
occupied by the chosen SNTF with a periodically nonhomogeneous permittivity 
dyadic [17, 18, 19] 

W^W = .oi^(.)-4,(z)-i;^W, (1) 

where the dyadics 

^^(z) = (u^Ux + UzVLz) cos [x{z)] + {nzVix - VixViz) sin [x{z)] + Uyiiy ^ 



gl^j:{z) = ea{z) UzUz + eb{z) u^u^ + e^z) UyU. 



(2) 



y"-v 



depend on the vapor incidence angle Xviz) = Xv + Sy sin(7rz/f2) that varies 
sinusoidally with period 2D,. 

In order to investigate SPP-wave propagation, we adopted a procedure de- 
vised to investigate the propagation of Dyakonov-Tamm waves [23]. Let the 
SPP wave propagate parallel to the unit vector u^; cos ip + Uy sin -0 along the 
interface z = 0, and attenuate as z — >■ ±oo. Therefore, in the region z < 0, the 
wave vector may be written as 

k„iet = K Ui - Umet ^z , (3) 

where + a^gj = k^ emet, ^ is complex- valued, and Im(amet) > for attenu- 
ation as z — — oo; here and hereafter, the unit vectors Ui = costp -f Uy simp 
and U2 = — sunp -\- Uy cosTp. Accordingly, the field phasors in the metal may 
be written as 



E(r) 



ui + — + as U2 



exp(ik„iet • r) , z <0, (4) 
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and 
H(r) = 77-1 



—o,p e„iet U2 + a. 



\ ^0 ^0 



exp(ikmet • r) , z <0, 



(5) 

where ap and are unknown scalars. 

For field representation in the SNTF, let us write E(r) = e(z) cxp (ikui • r) 
and H(r) = h(z) exp (ikui -r). The components e^iz) and hz{z) of the field 
phasors can be found in terms of the other components as follows: 



ed{z) [ea{z) - eb{z)] sin [x(z)] cos [xiz)] 



ea{z)ebiz) 



ex{z) 



^d{z) 



u}eoea{z)eb{z) 



[hx{z) sinijj — hy{z) cosip] , z>0. 



hz{z) 



[e2;(z) sin — ej,(z) 008-0] , z>0, 



where 



ea{z) eb{z) 



(6) 
(7) 

(8) 



eaiz) cos2 [xiz)] + ebiz) sin^ [xiz)] ' 
The other components of the electric and magnetic field phasors are used in 



the column vector 

[f(^)] = l^xiz) Cyiz) hxiz) hyiz)]^ 

which satisfies the matrix differential equation [23] 

^ [f (z)] = z ^(2)] • [f (z)] , z>0, 
where the 4x4 matrix 



[Piz)] 



(9) 
(10) 
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(11) 
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As in Part II, we used the piecewise uniform approximation technique [22] 
to determine the matrix [Q] that appears in the relation 



[f(2f])] = [Q] ■ [f(0+)] 



(12) 



to characterize the optical response of one period of the chosen SNTF for specific 
values of k and ip- By virtue of the Floquet-Lyapunov theorem [24], a matrix 
[Q] can be defined such that 



[£] =cxp{i2r![£]} 



(13) 



Both [Q] and [Q] share the same eigenvectors, and their eigenvalues are also 

related. Let [t]^"\ n e [1,4], be the eigenvector corresponding to the nth 
eigenvalue tT„ of [Q]; then, the corresponding eigenvalue a„ of [Q] is given by 



.In cr„ 



e[i,4] 



After ensuring that Im(Q!i_2) > 0, we set 

[f(o+)]- f[t](^) [t](^) 



b2 



(14) 



(15) 



for SPP-wave propagation, where hi and 62 a-re unknown scalars; the other two 
eigenvalues of [Q] pertain to waves that amplify as z — >■ cxd and cannot therefore 
contribute to the SPP wave. At the same time, 



[f(0- 



cos lb 
simh 

fen ^ 



tnijiL sin?/' 



vo 



— sin tJj 
cos-0 

cos^ 
- sinip 



(16) 



by virtue of (4) and (5). Continuity of the tangential components of the electric 
and magnetic field phasors across the plane z = requires that [f(0— )] = 
[f(0+)], which may be rearranged as the matrix equation 



Up 




" ' 









bi 















(17) 



For a nontrivial solution, the 4x4 matrix [Y] must be singular, so that 

det [r] = (18) 

is the dispersion equation for the SPP wave. This equation has to be solved in 
order to determine the SPP wavenumber k. 
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3 Numerical Results and Discussion 



A Mathcmatica'^'^ program was written and implemented to solve (18) using 
the Newton-Raphson method [26] to obtain k for a specific value of ip. The 
free-space wavelength was fixed at Ao = 633 nm. The metal was taken to be 
aluminum: tmet = ~56 + 21i. The SNTF was chosen to be made of titanium 
oxide [25], with 

ea{z) = [1.0443 + 2.7394t;(z) - 1.3697i;2(z)]2 ^ 
ebiz) = [1.6765 + 1.5649w(z) - Q.lS2bv^{z)Y 

) (19) 

ec(z) = [1.3586 + 2. 1109w(z) - 1.0554w2(z)]2 

X{z) tan-i[2.8818tanx„(2:)] 

where v{z) = 2xv{z)/tt. The angles Xv and 5y were taken to be 45° and 30°, 
respectively, and Vl = 200 nm for all results presented here. 
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Figure 1: (left) Real and (right) imaginary parts of k as functions of V', foi' SPP- 
wave propagation guided by the planar interface of aluminum and a titanium- 
oxide SNTF. Sec the text for the constitutive parameters used. Either two or 
three modes are possible, depending on tp. 
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The wavenumber k must be complex- valued for the canonical problem but 
only real- valued for the Kretschmann configuration [9] . Computed values of the 
real and imaginary parts of k for the canonical problem are shown in Fig. 1. 
For 0° < ^ 36°, we found three values of k which satisfy (18) and therefore 
represent SPP waves. For 36° ^ip < 90°, there are two values of k which satisfy 
(18). This trend is fully consistent with the conclusions drawn in Part II for 
the Kretschmann configuration. The different solutions of (18) for any specific 
value of ip indicate that the SPP waves have different phase speeds w/Re(K) — 
as theoretically predicted in Parts I and II, and experimentally confirmed in 
Part III — and different e-folding distances l/Im(K) along the direction of prop- 
agation. 

Specifically, for -0 = 0° the values of k which satisfy (18) are ki = (2.455 -I- 
i0.04208)fco, K2 = (2.080-l-i0.003538)fco, and Kg = (1.868 i0.007267)fco. These 
solutions represent SPP waves with wave vectors lying wholly in the morpho- 
logically significant plane of the SNTF, as addressed theoretically in Part I and 
experimentally in Part III. 
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Figure 2: Variations of components of e, h, and P with z, for k = (2.455 + 
i0.04208)A:o and ^ = 0°. The x-, y-, and 2-directed components are represented 
by sohd, dashed and chain-dashed hnes, respectively. The data were computed 
by setting Op = 1 V m~^, with Oj,, &i, and 62 then obtained using (17). 



The Cartesian components of the electric and magnetic field phasors and the 
time-averaged Poynting vector P = iRe (E x H*) as functions of z along the 
line {x = 0, y = 0) are shown for ki in Fig. 2, and for K3 in Fig. 4. These SPP 
waves are p-polarized, and were respectively labeled as pi and p2 in Part III. 
Figure 3 shows the variations of e, h, and P along the z axis for k = K2. This 
SPP wave is s-polarized and was labeled as s3 in Part III. The localization 
of all three SPP waves around the interface z = is evident from Figs. 2-4. 
Also, the SPP wave pi is more localized inside the SNTP than either p2 or s3. 
Furthermore, the phase speed of pi is higher than that of s3, which exceeds the 
phase speed of p2. However, s3 will travel a longer distance along the interface 
than either pi or p2, which could not have been deduced from the theoretical 
analysis for the Kretschmann configuration in Part I. 

For -0 = 75°, only two values of k were found to satisfy (18): ki = (2.459 + 
i0.04247)fco and K2 = (2.066 + i0.003861)A:o. The variations of e, h, and P 
along the z axis are shown in Fig. 5 for ki, and in Fig. 6 for K2. These two 
SPP waves cannot be classified as either p- or s-polarized, which is consistent 
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Figure 3: Same as Fig. 2 except for k = (2.080 + i0.003538)A:o. The data were 
computed by setting Os = 1 V m~^, with Us, bi, and 62 then obtained using 
(17). 



with the deductions in Part II for the Kretschmann configuration. Although 
the phase speed of the SPP wave with k — ki exceeds the phase speed of the 
other SPP wave {k = K2), the latter SPP wave will propagate a longer distance 
along the interface than the former SPP wave. 

4 Concluding Remarks 

In Parts I and II of this series, the possibility of excitation of multiple SPP waves 
guided by the planar interface of a metal and a periodically nonhomogeneous 
SNTP was predicted after solving a boundary-value problem that models the 
Kretschmann setup. In Part III, an experimental validation was provided for 
the results of Part I. However, these approaches were indirect and a direct proof 
of the existence of multiple SPP waves was lacking. 

The direct proof was provided here after formulating and solving the canon- 
ical problem of propagation guided by the planar interface of semi-infinite ex- 
panses of the metal and the periodically nonhomogeneous SNTP. Not only were 
the conclusions obtained in Parts I-III upheld, but additional information on 
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the e-folding distance along the direction of propagation also emerged from the 
solution of the canonical problem. 
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